Report for Model 5

— for both jointly and separately

2024-08-10

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.02
## LDevsum       1.01       1.03
## dh0           1.50       2.15
## dh1           1.66       2.63
## dh2           1.41       2.22
## dl0           1.92       3.06
## dl1           1.60       2.49
## dl2           1.39       2.11
## 
## Multivariate psrf
## 
## 1.83

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1158.7199 22017.5632
DIC3 1122.8126 22085.9396
PWAIC 26.4105 187.0915
WAIC 1136.2745 22098.6314

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.00       1.01
## LDevsum       1.00       1.01
## dh0           2.55       8.40
## dh1           9.41      17.85
## dh2           5.53      13.34
## dl0           1.12       1.33
## dl1           1.17       1.47
## dl2           1.05       1.14
## 
## Multivariate psrf
## 
## 8.43

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1380.41175 30292.6188
DIC3 1244.84499 30382.9877
PWAIC 90.14513 139.2839
WAIC 1313.83189 30387.0876